Differential gene expression in layer 3 pyramidal neurons across 3 regions of the human cortical visual spatial working memory network.

Visual spatial working memory (vsWM) is mediated by a distributed cortical network composed of multiple nodes, including primary visual (V1), posterior parietal (PPC), and dorsolateral prefrontal (DLPFC) cortices. Feedforward and feedback information is transferred among these nodes via projections furnished by pyramidal neurons (PNs) located primarily in cortical layer 3. Morphological and electrophysiological differences among layer 3 PNs across these nodes have been reported; however, the transcriptional signatures underlying these differences have not been examined in the human brain. Here we interrogated the transcriptomes of layer 3 PNs from 39 neurotypical human subjects across 3 critical nodes of the vsWM network. Over 8,000 differentially expressed genes were detected, with more than 6,000 transcriptional differences present between layer 3 PNs in V1 and those in PPC and DLPFC. Additionally, over 600 other genes differed in expression along the rostral-to-caudal hierarchy formed by these 3 nodes. Moreover, pathway analysis revealed enrichment of genes in V1 related to circadian rhythms and in DLPFC of genes involved in synaptic plasticity. Overall, these results show robust regional differences in the transcriptome of layer 3 PNs, which likely contribute to regional specialization in their morphological and physiological features and thus in their functional contributions to vsWM.

Lung transplantation for idiopathic pulmonary fibrosis enriches for individuals with telomere-mediated disease.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is the most common indication for lung transplantation in North America and variants in telomere-maintenance genes are the most common identifiable cause of IPF. We reasoned that younger IPF patients are more likely to undergo lung transplantation and we hypothesized that lung transplant recipients would be enriched for individuals with telomere-mediated disease due to the earlier onset and more severe disease in these patients. METHODS: Individuals with IPF who underwent lung transplantation or were evaluated in an interstitial lung disease specialty clinic who did not undergo lung transplantation were examined. Genetic evaluation was completed via whole genome sequencing (WGS) of 426 individuals and targeted sequencing for 5 individuals. Rare variants in genes previously associated with IPF were classified using the American College of Medical Genetics guidelines. Telomere length from WGS data was measured using TelSeq software. Patient characteristics were collected via medical record review. RESULTS: Of 431 individuals, 149 underwent lung transplantation for IPF. The median age of diagnosis of transplanted vs non-transplanted individuals was significantly younger (60 years vs 70 years, respectively, p<0.0001). IPF lung transplant recipients (IPF-LTRs) were twice as likely to have telomere-related rare variants compared to non-transplanted individuals (24% vs 12%, respectively, p=0.0013). IPF-LTRs had shorter telomeres than non-transplanted IPF patients (p=0.0028) and >85% had telomeres below the age-adjusted mean. Post-transplant survival and CLAD were similar amongst IPF-LTRs with rare variants in telomere-maintenance genes compared to those without, as well as in those with short telomeres versus longer telomeres. CONCLUSIONS: There is an enrichment for telomere-maintenance gene variants and short telomeres among IPF-LTRs. However, transplant outcomes of survival and CLAD do not differ by gene variants or telomere length within IPF-LTRs. Our findings support individual with telomere-mediated disease should not be excluded from lung transplantation and focusing research efforts on therapies directed toward individuals with short-telomere mediated disease.

A Genome-Wide Analysis of the Penumbral Volume in Inbred Mice following Middle Cerebral Artery Occlusion.

Following ischemic stroke, the penumbra, at-risk neural tissue surrounding the core infarct, survives for a variable period of time before progressing to infarction. We investigated genetic determinants of the size of penumbra in mice subjected to middle cerebral artery occlusion (MCAO) using a genome-wide approach. 449 male mice from 33 inbred strains underwent MCAO for 6 hours (215 mice) or 24 hours (234 mice). A genome-wide association study using genetic data from the Mouse HapMap project was performed to examine the effects of genetic variants on the penumbra ratio, defined as the ratio of the infarct volume after 6 hours to the infarct volume after 24 hours of MCAO. Efficient mixed model analysis was used to account for strain interrelatedness. Penumbra ratio differed significantly by strain (F = 2.7, P < 0.001) and was associated with 18 significant SNPs, including 6 protein coding genes. We have identified 6 candidate genes for penumbra ratio: Clint1, Nbea, Smtnl2, Rin3, Dclk1, and Slc24a4.

Variable Susceptibility to Cigarette Smoke-Induced Emphysema in 34 Inbred Strains of Mice Implicates Abi3bp in Emphysema Susceptibility.

Chronic obstructive pulmonary disease (COPD) is caused by a complex interaction of environmental exposures, most commonly cigarette smoke, and genetic factors. Chronic cigarette smoke exposure in the mouse is a commonly used animal model of COPD. We aimed to expand our knowledge about the variable susceptibility of inbred strains to this model and test for genetic variants associated with this trait. To that end, we sought to measure differential susceptibility to cigarette smoke-induced emphysema in the mouse, identify genetic loci associated with this quantitative trait, and find homologous human genes associated with COPD. Alveolar chord length (CL) in 34 inbred strains of mice was measured after 6 months of exposure to cigarette smoke. After testing for association, we connected a murine candidate locus to a published meta-analysis of moderate-to-severe COPD. We identified deleterious mutations in a candidate gene in silico and measured gene expression in extreme strains. A/J was the most susceptible strain in our survey (Δ CL 7.0 ± 2.2 µm) and CBA/J was the least susceptible (Δ CL -0.3 ± 1.2 µm). By integrating mouse and human genome-wide scans, we identified the candidate gene Abi3bp. CBA/J mice harbor predicted deleterious variants in Abi3bp, and expression of the gene differs significantly between CBA/J and A/J mice. This is the first report of susceptibility to cigarette smoke-induced emphysema in 34 inbred strains of mice, and Abi3bp is identified as a potential contributor to this phenotype.

Genome-Wide Association Mapping of the Antibody Response to Diphtheria, Tetanus and Acellular Pertussis Vaccine in Mice.

The objective of the current study was to investigate the genetics of antibody responses to an acellular pertussis vaccine by a genome-wide association study in mice. Female mice of 28 inbred strains received this vaccine at 6, 8, and 12 weeks of age. The antibody titer and avidity of immunoglobulin (Ig) G specific for diphtheria toxin, pertussis toxin, filamentous hemagglutinin and pertactin were measured at 14 and 24 weeks of age. The magnitude, longevity and avidity of IgG differed significantly among mouse strains. There was significant correlation between antigen-specific IgGs for longevity but not for magnitude and avidity. Association mapping and analysis with PolyPhen software identified 6 genetic markers associated with longevity for all 4 antigens, although the expression levels of these genes did not correlate with longevity phenotype. This study provides novel insights into the genetic basis and potential candidate genes for differences in the IgG responses to vaccination.

TCGA Expedition: A Data Acquisition and Management System for TCGA Data.

BACKGROUND: The Cancer Genome Atlas Project (TCGA) is a National Cancer Institute effort to profile at least 500 cases of 20 different tumor types using genomic platforms and to make these data, both raw and processed, available to all researchers. TCGA data are currently over 1.2 Petabyte in size and include whole genome sequence (WGS), whole exome sequence, methylation, RNA expression, proteomic, and clinical datasets. Publicly accessible TCGA data are released through public portals, but many challenges exist in navigating and using data obtained from these sites. We developed TCGA Expedition to support the research community focused on computational methods for cancer research. Data obtained, versioned, and archived using TCGA Expedition supports command line access at high-performance computing facilities as well as some functionality with third party tools. For a subset of TCGA data collected at University of Pittsburgh, we also re-associate TCGA data with de-identified data from the electronic health records. Here we describe the software as well as the architecture of our repository, methods for loading of TCGA data to multiple platforms, and security and regulatory controls that conform to federal best practices. RESULTS: TCGA Expedition software consists of a set of scripts written in Bash, Python and Java that download, extract, harmonize, version and store all TCGA data and metadata. The software generates a versioned, participant- and sample-centered, local TCGA data directory with metadata structures that directly reference the local data files as well as the original data files. The software supports flexible searches of the data via a web portal, user-centric data tracking tools, and data provenance tools. Using this software, we created a collaborative repository, the Pittsburgh Genome Resource Repository (PGRR) that enabled investigators at our institution to work with all TCGA data formats, and to interrogate these data with analysis pipelines, and associated tools. WGS data are especially challenging for individual investigators to use, due to issues with downloading, storage, and processing; having locally accessible WGS BAM files has proven invaluable. CONCLUSION: Our open-source, freely available TCGA Expedition software can be used to create a local collaborative infrastructure for acquiring, managing, and analyzing TCGA data and other large public datasets.

Mouse genome-wide association study identifies polymorphisms on chromosomes 4, 11, and 15 for age-related cardiac fibrosis.

Dystrophic cardiac calcinosis (DCC), also called epicardial and myocardial fibrosis and mineralization, has been detected in mice of a number of laboratory inbred strains, most commonly C3H/HeJ and DBA/2J. In previous mouse breeding studies between these DCC susceptible and the DCC-resistant strain C57BL/6J, 4 genetic loci harboring genes involved in DCC inheritance were identified and subsequently termed Dyscalc loci 1 through 4. Here, we report susceptibility to cardiac fibrosis, a sub-phenotype of DCC, at 12 and 20 months of age and close to natural death in a survey of 28 inbred mouse strains. Eight strains showed cardiac fibrosis with highest frequency and severity in the moribund mice. Using genotype and phenotype information of the 28 investigated strains, we performed genome-wide association studies (GWAS) and identified the most significant associations on chromosome (Chr) 15 at 72 million base pairs (Mb) (P < 10(-13)) and Chr 4 at 122 Mb (P < 10(-11)) and 134 Mb (P < 10(-7)). At the Chr 15 locus, Col22a1 and Kcnk9 were identified. Both have been reported to be morphologically and functionally important in the heart muscle. The strongest Chr 4 associations were located approximately 6 Mb away from the Dyscalc 2 quantitative trait locus peak within the boundaries of the Extl1 gene and in close proximity to the Trim63 and Cap1 genes. In addition, a single-nucleotide polymorphism association was found on chromosome 11. This study provides evidence for more than the previously reported 4 genetic loci determining cardiac fibrosis and DCC. The study also highlights the power of GWAS in the mouse for dissecting complex genetic traits.

Genetic determinants of fibro-osseous lesions in aged inbred mice.

Fibro-osseous lesions in mice are progressive aging changes in which the bone marrow is replaced to various degrees by fibrovascular stroma and bony trabeculae in a wide variety of bones. The frequency and severity varied greatly among 28 different inbred mouse stains, predominantly affecting females, ranging from 0% for 10 strains to 100% for KK/HlJ and NZW/LacJ female mice. Few lesions were observed in male mice and for 23 of the strains, no lesions were observed in males for any of the cohorts. There were no significant correlations between strain-specific severities of fibro-osseous lesions and ovarian (r=0.11; P=0.57) or endometrial (r=0.03; P=0.89) cyst formation frequency or abnormalities in parathyroid glands. Frequency of fibro-osseous lesions was most strongly associated (P<10(-6)) with genome variations on chromosome (Chr) 8 at 90.6 and 90.8Mb (rs33108071, rs33500669; P=5.0·10(-10), 1.3·10(-6)), Chr 15 at 23.6 and 23.8Mb (rs32087871, rs45770368; P=7.3·10(-7), 2.7·10(-6)), and Chr 19 at 33.2, 33.4, and 33.6Mb (rs311004232, rs30524929, rs30448815; P=2.8·10(-6), 2.8·10(-6), 2.8·10(-6)) in genome-wide association studies (GWAS). The relatively large number of candidate genes identified in the GWAS analyses suggests that this may be an extremely complex polygenic disease. These results indicate that fibro-osseous lesions are surprisingly common in many inbred strains of laboratory mice as they age. While this presents little problem in most studies that utilize young animals, it may complicate aging studies, particularly those focused on bone.

Integrative Mouse and Human Studies Implicate ANGPT1 and ZBTB7C as Susceptibility Genes to Ischemic Injury.

BACKGROUND AND PURPOSE: The extent of ischemic injury in response to cerebral ischemia is known to be affected by native vasculature. However, the nonvascular and dynamic vascular responses and their genetic basis are not well understood. METHODS: We performed a genome-wide association study in 235 mice from 33 inbred strains using the middle cerebral artery occlusion model. Population structure and genetic relatedness were accounted for using the efficient mixed-model association method. Human orthologs to the genes associated with the significant and suggestive single-nucleotide polymorphisms from the mouse strain survey were examined in patients with M1 occlusions admitted with signs and symptoms of acute ischemic stroke. RESULTS: We identified 4 genome-wide significant and suggestive single-nucleotide polymorphisms to be associated with infarct volume in mice (rs3694965, P=2.17×10(-7); rs31924033, P=5.61×10(-6); rs32249495, P=2.08×10(-7); and rs3677406, P=9.56×10(-6)). rs32249495, which corresponds to angiopoietin-1 (ANGPT1), was also significant in the recessive model in humans, whereas rs1944577, which corresponds to ZBTB7C, was nominally significant in both the additive and dominant genetic models in humans. ZBTB7C was shown to be upregulated in endothelial cells using both in vitro and in vivo models of ischemia. CONCLUSIONS: Genetic variations of ANGPT1 and ZBTB7C are associated with increased infarct size in both mice and humans. ZBTB7C may modulate the ischemic response via neuronal apoptosis and dynamic collateralization and, in addition to ANGPT1, may serve as potential novel targets for treatments of cerebral ischemia.

Mouse models for pseudoxanthoma elasticum: genetic and dietary modulation of the ectopic mineralization phenotypes.

Pseudoxanthoma elasticum (PXE), a heritable ectopic mineralization disorder, is caused by mutations in the ABCC6 gene. Null mice (Abcc6(-/-) ) recapitulate the genetic, histopathologic and ultrastructural features of PXE, and they demonstrate early and progressive mineralization of vibrissae dermal sheath, which serves as a biomarker of the overall mineralization process. Recently, as part of a mouse aging study at The Jackson Laboratory, 31 inbred mouse strains were necropsied, and two of them, KK/HlJ and 129S1/SvImJ, were noted to have vibrissae dermal mineralization similar to Abcc6(-/-) mice. These two strains were shown to harbor a single nucleotide polymorphism (rs32756904) in the Abcc6 gene, which resulted in out-of-frame splicing and marked reduction in ABCC6 protein expression in the liver of these mice. The same polymorphism is present in two additional mouse strains, DBA/2J and C3H/HeJ, with similar reduction in Abcc6 protein levels, yet these mice did not demonstrate tissue mineralization when kept on standard rodent diet. However, all four mouse strains, when placed on experimental diet enriched in phosphate and low in magnesium, developed extensive ectopic mineralization. These results indicate that the genetic background of mice and the mineral composition of their diet can profoundly modulate the ectopic mineralization process predicated on mutations in the Abcc6 gene. These mice provide novel model systems to study the pathomechanisms and the reasons for strain background on phenotypic variability of PXE.

Personalized medicine for chronic, complex diseases: chronic obstructive pulmonary disease as an example.

Chronic, complex diseases represent the majority of healthcare utilization and spending in the USA today. Despite this, therapeutics that account for the heterogeneity of these diseases are lacking, begging for more personalized approaches. Improving our understanding of disease phenotypes through retrospective trials of electronic health record data will enable us to better categorize patients. Increased usage of next-generation sequencing will further our understanding of the genetic variants involved in chronic disease. Utilization of data warehousing will be necessary in order to securely handle, integrate and analyze the large sets of data produced with these methods. Finally, increased use of clinical decision support will enable the return of clinically actionable results that physicians can use to apply these personalized approaches.

Alopecia areata: updates from the mouse perspective.

Alopecia areata (AA) is a cell-mediated autoimmune disease that targets actively growing hair follicles in mammals, including humans and mice. Development of the C3H/HeJ spontaneous mouse model AA nearly 20 years ago provided a much needed tool to test the hypotheses and ultimately serve as a preclinical model for drug testing. Discoveries in both human AA patients and the mouse model supported each other and lead to discoveries on the incredibly complex genetic basis of this disease. The discovery that A/J, MRL/MpJ, SJL/J, and SWR/J strains also develop AA now allows genome-wide association mapping studies to expand the list of genes underlying this disease. Potential new targets for unraveling the pathogenesis of AA include the role of retinoic acid metabolism in the severity of disease and hair shaft proteins that may be either the inciting antigen or ultimate target of the immune reaction leading to breakage of the shaft causing clinical alopecia. Comparing these model systems with human and mouse clinical disease, for both discovery and validation of the discoveries, continues to resolve the complex questions surrounding AA.

Mutant Enpp1asj mice as a model for generalized arterial calcification of infancy.

Generalized arterial calcification of infancy (GACI), an autosomal recessive disorder, is characterized by early mineralization of blood vessels, often diagnosed by prenatal ultrasound and usually resulting in demise during the first year of life. It is caused in most cases by mutations in the ENPP1 gene, encoding an enzyme that hydrolyzes ATP to AMP and inorganic pyrophosphate, the latter being a powerful anti-mineralization factor. Recently, a novel mouse phenotype was recognized as a result of ENU mutagenesis - those mice developed stiffening of the joints, hence the mutant mouse was named 'ages with stiffened joints' (asj). These mice harbor a missense mutation, p.V246D, in the Enpp1 gene. Here we demonstrate that the mutant ENPP1 protein is largely absent in the liver of asj mice, and the lack of enzymatic activity results in reduced inorganic pyrophosphate (PPi) levels in the plasma, accompanied by extensive mineralization of a number of tissues, including arterial blood vessels. The progress of mineralization is highly dependent on the mineral composition of the diet, with significant shortening of the lifespan on a diet enriched in phosphorus and low in magnesium. These results suggest that the asj mouse can serve as an animal model for GACI.

Functional genomic assessment of phosgene-induced acute lung injury in mice.

In this study, a genetically diverse panel of 43 mouse strains was exposed to phosgene and genome-wide association mapping performed using a high-density single nucleotide polymorphism (SNP) assembly. Transcriptomic analysis was also used to improve the genetic resolution in the identification of genetic determinants of phosgene-induced acute lung injury (ALI). We prioritized the identified genes based on whether the encoded protein was previously associated with lung injury or contained a nonsynonymous SNP within a functional domain. Candidates were selected that contained a promoter SNP that could alter a putative transcription factor binding site and had variable expression by transcriptomic analyses. The latter two criteria also required that ≥10% of mice carried the minor allele and that this allele could account for ≥10% of the phenotypic difference noted between the strains at the phenotypic extremes. This integrative, functional approach revealed 14 candidate genes that included Atp1a1, Alox5, Galnt11, Hrh1, Mbd4, Phactr2, Plxnd1, Ptprt, Reln, and Zfand4, which had significant SNP associations, and Itga9, Man1a2, Mapk14, and Vwf, which had suggestive SNP associations. Of the genes with significant SNP associations, Atp1a1, Alox5, Plxnd1, Ptprt, and Zfand4 could be associated with ALI in several ways. Using a competitive electrophoretic mobility shift analysis, Atp1a1 promoter (rs215053185) oligonucleotide containing the minor G allele formed a major distinct faster-migrating complex. In addition, a gene with a suggestive SNP association, Itga9, is linked to transforming growth factor ß1 signaling, which previously has been associated with the susceptibility to ALI in mice.

Integration of mouse and human genome-wide association data identifies KCNIP4 as an asthma gene.

Asthma is a common chronic respiratory disease characterized by airway hyperresponsiveness (AHR). The genetics of asthma have been widely studied in mouse and human, and homologous genomic regions have been associated with mouse AHR and human asthma-related phenotypes. Our goal was to identify asthma-related genes by integrating AHR associations in mouse with human genome-wide association study (GWAS) data. We used Efficient Mixed Model Association (EMMA) analysis to conduct a GWAS of baseline AHR measures from males and females of 31 mouse strains. Genes near or containing SNPs with EMMA p-values <0.001 were selected for further study in human GWAS. The results of the previously reported EVE consortium asthma GWAS meta-analysis consisting of 12,958 diverse North American subjects from 9 study centers were used to select a subset of homologous genes with evidence of association with asthma in humans. Following validation attempts in three human asthma GWAS (i.e., Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG) and two human AHR GWAS (i.e., SHARP, DAG), the Kv channel interacting protein 4 (KCNIP4) gene was identified as nominally associated with both asthma and AHR at a gene- and SNP-level. In EVE, the smallest KCNIP4 association was at rs6833065 (P-value 2.9e-04), while the strongest associations for Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG were 1.5e-03, 1.0e-03, 3.1e-03 at rs7664617, rs4697177, rs4696975, respectively. At a SNP level, the strongest association across all asthma GWAS was at rs4697177 (P-value 1.1e-04). The smallest P-values for association with AHR were 2.3e-03 at rs11947661 in SHARP and 2.1e-03 at rs402802 in DAG. Functional studies are required to validate the potential involvement of KCNIP4 in modulating asthma susceptibility and/or AHR. Our results suggest that a useful approach to identify genes associated with human asthma is to leverage mouse AHR association data.

Early gene expression differences in inbred mouse strains with susceptibility to pulmonary adenomas.

Lung cancer is the most common cause of cancer-related deaths in both men and women, and effective preventatives are rare due to the difficulty of early detection. Specific gene expression signatures have been identified in individuals that already developed lung cancer. To identify if gene expression differences could be detected in individuals before the onset of the disease, we obtained lung tissues for microarray analysis from young, healthy mice of 9 inbred strains with known differences in their susceptibility to spontaneous pulmonary adenomas when aged. We found that the most common differentially expressed genes among all possible 36 strain comparisons showed significant associations with cancer- and inflammation-related processes. Significant expression differences between susceptible and resistant strains were detected for Aldh3a1, Cxcr1 and 7, Dpt, and Nptx1-genes with known cancer-related functions, and Cd209, Cxcr1 and 7, and Plag2g1b-genes with known inflammatory-related functions. Whereas Aldh3a1, Cd209, Dpt, and Pla2g1b had increased expression, Cxcr1 and 7, and Nptx1 had decreased expression in strains susceptible to pulmonary adenomas. Thus, our study shows that expression differences between susceptible and resistant strains can be detected in young and healthy mice without manifestation of pulmonary adenomas and, thus, may provide an opportunity of early detection. Finally, the identified genes have previously been reported for human non-small cell lung cancer suggesting that molecular pathways may be shared between these two cancer types.

Emerging genetics of COPD.

Since the discovery of alpha-1 antitrypsin in the early 1960s, several new genes have been suggested to play a role in chronic obstructive pulmonary disease (COPD) pathogenesis. Yet, in spite of those advances, much about the genetic basis of COPD still remains to be discovered. Unbiased approaches, such as genome-wide association (GWA) studies, are critical to identify genes and pathways and to verify suggested genetic variants. Indeed, most of our current understanding about COPD candidate genes originates from GWA studies. Experiments in form of cross-study replications and advanced meta-analyses have propelled the field towards unravelling details about COPD's pathogenesis. Here, we review the discovery of genetic variants in association with COPD phenotypes by discussing the available approaches and current findings. Limitations of current studies are considered and future directions provided.

A major X-linked locus affects kidney function in mice.

Chronic kidney disease is a common disease with increasing prevalence in the western population. One common reason for chronic kidney failure is diabetic nephropathy. Diabetic nephropathy and hyperglycemia are characteristics of the mouse inbred strain KK/HlJ, which is predominantly used as a model for metabolic syndrome due to its inherited glucose intolerance and insulin resistance. We used KK/HlJ, an albuminuria-sensitive strain, and C57BL/6J, an albuminuria-resistant strain, to perform a quantitative trait locus (QTL) cross to identify the genetic basis for chronic kidney failure. Albumin-creatinine ratio (ACR) was measured in 130 F2 male offspring. One significant QTL was identified on chromosome (Chr) X and four suggestive QTL were found on Chrs 6, 7, 12, and 13. Narrowing of the QTL region was focused on the X-linked QTL and performed by incorporating genotype and expression analyses for genes located in the region. From the 485 genes identified in the X-linked QTL region, a few candidate genes were identified using a combination of bioinformatic evidence based on genomic comparison of the parental strains and known function in urine homeostasis. Finally, this study demonstrates the significance of the X chromosome in the genetic determination of albuminuria.